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Abstract 

Cross sections for semi-inclusive electroproduction of charged pions (tt^) from both proton and 
deuteron targets were measured for 0.2 < x < 0.5, 2 < Q 2 < 4 GeV 2 , 0.3 < z < 1, and Pj 2 < 0.2 
GeV 2 . For P t < 0.1 GeV, we find the azimuthal dependence to be small, as expected theoretically. 
For both ir + and 7r _ , the Pt dependence from the deuteron is found to be slightly weaker than from 
the proton. In the context of a simple model, this implies that the initial transverse momenta width 
of d quarks is larger than for u quarks and, contrary to expectations, the transverse momentum 
width of the favored fragmentation function is larger than the unfavored one. 

PACS numbers: 13.60.Le, 13.87.Fn 
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FIG. 1: Schematic diagram of semi- inclusive pion electroproduction within a factorized QCD 
parton model at lowest order in a s . Final transverse momenta of the detected pion Pt arises 
from convolving the struck quark transverse momenta kt with the transverse momentum generated 
during fragmentation process pt- 



A central question in the understanding of nucleon structure is the orbital motion of 
partons. Much is known about the light-cone momentum fraction, x, and virtuality scale, 
Q 2 , dependence of the up and down quark parton distribution functions (PDFs) in the 
nucleon. In contrast, very little is presently known about the dependence of these functions 
on their transverse momentum k t . Simply based on the size of the nucleon in which the 
quarks are confined, one would expect characteristic transverse momenta of order a few 
hundred MeV, with larger values at small Bjorken x where the sea quarks dominate, and 
smaller values at high x where all of the quark momentum is longitudinal in the limit x — 1. 
Increasingly precise studies of the nucleon spin sum rule , |4| strongly suggest that the 

net spin carried by quarks and gluons is relatively small, and therefore the net orbital angular 
momentum must be significant. This in turn implies significant transverse momentum of 
quarks. Questions that naturally arise include: what is the flavor and helicity dependence 
of the transverse motion of quarks and gluons, and can these be modeled theoretically and 
measured experimentally? 

The process of semi- inclusive deep- inelastic lepton scattering (SIDIS), IN — > IhX has 
been shown to f actor ize jsj], in the high energy limit, into lepton-quark scattering followed 
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by quark hadronization. Ideally, one could directly measure the quark transverse momentum 
dependence of the quark distribution functions q(x, k t ) by detecting all particles produced in 
the hadronization process. In the present experiment, we detect only a single hadronization 
product: a charged pion carrying an energy fraction z of the available energy. The probability 
of producing a pion with a transverse momentum P t relative to the virtual photon (q) 
direction is described by a convolution of the quark distribution functions and ^-dependent 
fragmentation functions D + (z,p t ) and D~(z,p t ), where pt is the transverse momentum of 
the pion relative to the quark direction, with the imposed condition [6] P t = zk t + p t 
(see Fig. [1]). The "favored" and "unfavored" functions D + (z,p t ) and D~(z,p t ) refer to 
the case where the produced pion contains of the same flavor as the struck quark or not. 
"Soft" non-perturbative processes are expected 6| to generate relatively small values of p t 
with an approximately Gaussian distributions in p t . Hard QCD processes are expected to 
generate large non-Gaussian tails for p t > 1 GeV, and probably do not play a major role 
in the interpretation of the present experiment, for which the total transverse momentum 
Pt < 0.45 GeV. The assumption that the fragmentation functions do not depend on quark 
flavor (for example D + (z,pt) applies equally well to u — > ti + and d — > 7r~) in principle allows 
the kt widths of up and down quarks to be distinguished. In the present experiment, the 
use of both proton and deuteron targets (the latter with a higher d quark content than the 
former) and the detection of both ir + and tt~ permits a first study of this problem. 

The experiment (E00-108) used the Short Orbit (SOS) and High Momentum (HMS) 
spectrometers in Hall C at Jefferson Lab to detect final state electrons and pions, respectively. 
An electron beam with energy of 5.5 GeV and currents ranging between 20 and 60 \xA was 
provided by the CEBAF accelerator. Incident electrons were scattered from 4-cm-long liquid 
hydrogen or deuterium targets. The experiment consisted of three parts: i) at a fixed electron 
kinematics of (x, Q 2 ) = (0.32, 2.30 GeV 2 ), z was varied from 0.3 to 1, with nearly uniform 
coverage in the pion azimuthal angle, 0, around the virtual photon direction, but at a small 
average Pt of 0.05 GeV; ii) for z = 0.55, x was varied from 0.2 to 0.5 (with a corresponding 
variation in Q 2 , from 2 to 4 GeV 2 ), keeping the pion centered on the virtual photon direction 
(and again average P t of 0.05 GeV); iii) for (x, Q 2 ) = (0.32, 2.30 GeV 2 ), z near 0.55, P t was 
scanned from to 0.4 GeV by increasing the HMS angle (with average near 180 degrees). 
The <p distribution as a function of P t is shown for all three data sets combined in Fig. [21 
The virtual photon-nucleon invariant mass W, was always larger than 2.1 GeV (typically 



4 



0.5 




50 100 150 200 250 300 350 



* (deg) 



FIG. 2: Pi distribution of data from this experiment as a function of 0. 

2.4 GeV), corresponding to the traditional deep inelastic region for inclusive scattering. 

At lower virtual photon energy and/or mass scales, the factorization ansatz is expected 
to break down, due to the effects of final state interactions, resonant nucleon excitations, 
and higher twist contributions [7]. In particular, in the present experiment the residual 
invariant mass M x of the undetected particles (see Fig. [T|) ranges from about 1 to 2 GeV 
(inversely correlated with z), spanning the mass region traditionally associated with signifi- 
cant baryon resonance excitation. The extent to which this situation leads to a break-down 
of factorization was studied in our previous paper jg]. It was found that good agreement 
with expectations based on higher energy data was achieved for z < 0.7, approximately cor- 
responding to M x > 1.5 GeV. The ratio of total up to down quark distributions u(x)/d(x) 
extracted from ratios of cross sections, as well as the ratio of valence-only up to down ratios 
u v (x)/d v (x), were also found to be reasonably compatible with higher energy extractions, 
provided z < 0.7. This issue will be addressed further for the Pf-scan data below. Finally, 
the ratio of unfavored to favored fragmentation functions D~(z)/D + (z) (from the 7r~/7r + 
ratios on the deuteron) was found to be consistent with extractions from other experiments. 
All of these studies were done with the z-scan and x-scan data, for which the average Pt 
was small (< 0.1 GeV ), and the average value of cos(0) was close to zero. 

In this paper, we focus on the P t dependence, with the goal of searching for a possible 
flavor dependence to the quark distribution functions and/or fragmentation functions. Since 
the average value of cos(0) in the present experiment is correlated with P t (approaching -1 
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for the largest Pt value of 0.45 GeV, see Fig. [2]), we first study the limited data available 
from this experiment on the <fi dependence, which must be an even function since neither 
the beam nor the target were polarized. We parameterize the data for each target and 
pion flavor according to: 

da ee ^ x dN 2 1 + A cos <f) + B cos(20) 

— ; — — i— b expl — bP f ) (1) 

where the parameters A(x, Q 2 , z, P t ) and B(x,Q 2 ,z,Pt) are a measure of the relative im- 
portance of the interference terms <tlt and cttt, respectively [h]]. The assumed Gaussian 
P 2 dependence (with slopes b for each case) is an effective parameterization that seems to 
describe the data adequately for use in making radiative and bin-centering corrections. We 
use this model for studying the <p dependence, then return to a more detailed study of the 
Pt dependence in the context of a simple model that incorporates a different Pt dependence 
for each struck quark and produced hadron flavor. 

For each kinematic point in the x and z scans (average P t = 0.05 GeV, maximum P t 0.2 
GeV), we extracted A and B and found no statistically significant difference between the 
results for 7r + or 7r~, or proton or deuteron targets. We therefore combined all four cases 
together, and present the results in Fig. [31 Systematic errors (not shown in the figure) are 
approximately 0.03 on both A and B and are highly correlated from point to point. Taking 
the systematic errors into account, the values of A and B are close to zero, for all values of 
x studied, and for values of z < 0.7, where our previous studies showed a good consistency 
with factorization. The small values of A and B are also consistent with the expectations 
based on kinematic shifts due to parton motion as described by Cahn [111] (shown as the solid 
curves on the figures) and Levelt-Mulders [3] . These effects are proportional to P t for A, and 
P 2 for B respectively [ll|, Q, [l3, JjJ], so are suppressed at low P t . More specifically, using 



the assumption that the average quark and fragmentation widths are equal, the Cahn 11 | 
asymmetries are given by 



A = - 7 (2 < Pt > /Q)(2 - y)y/T=y/[l + (1 - y)\ (2) 
B = - 7 2 (2 < P 2 t > /Q 2 )(l - y)/[l + (1 - y)\ (3) 

where 7 = z 2 / (1 + z 2 ), y — v/E, v is the virtual photon energy, and E is the beam energy, 
yielding A = —0.01 and B = —0.0002 for z = 0.55. The more recent treatment of Ref. [6| 
also gives results for A and B which are very close to zero (especially for B) . Other possible 



higher twist contributions will also be proportional to powers of Pt/^ftQ 2 ) m m, and 
therefore suppressed at our lower average values of P t and P 2 . Specifically, the twist- 2 
Boer-Mulders [l^] contribution to B is essentially zero in the models of Ref. 



17 



18|- 



In contrast, the longitudinal-transverse and transverse-transverse coefficients A and B 
are much larger in exclusive pion production (M x = M, where M is the nucleon mass) 
than those predicted for SIDIS. This is evidenced by our extracted average values for ex- 
clusive 7r electroproduction on deuteron and for ir + on proton, shown as the open symbols 
near z = 0.98 in Fig. [3j This underlines the importance of accounting for the radiative tail 
from exclusive production, which in our analysis was done using the computer code EX- 
CLURAD Q together with a reasonable model of exclusive pion electroproduction. The 
corrections where checked with the Hall C simulation package SIMC, which treats radiative 
corrections in the energy and angle peaking approximation 20j. 

We now turn to the study of the P t scan data. We used the cross section model from 
our previous paper [8j to describe the Q 2 dependence of the data (needed because P t and 
Q 2 are somewhat correlated), and extracted cross sections at fixed Q 2 averaged over 0. 
The corrections for Q 2 dependence did not distinguish between targets or pion flavor. Rel- 
atively small corrections (typically a few percent) for radiative effects (including the tails 
from exclusive pion production) and diffractive p production were made for each case 
individually. The systematic error on these corrections is estimated to be approximately 2%. 
The normalization errors due to target thickness, computer and electronic dead time, beam 
charge measurement, beam energy, and spectrometer kinematics combine to approximately 
2% overall, and 1% from case to case. The overall error due to spectrometer acceptance is 
estimated to be 3%, but < 1% from case to case because targets were exchanged frequently, 
as was the spectrometer polarity. The extracted cross sections are shown in Fig. H] and listed 
in Table [H The acceptance-averaged values of cos(0) range from -0.3 at low P t to nearly -1 
at higher P t , while the average values of cos(20) approaches 1 at high P t (See Fig. [2] and 
Table!]). 

Examination of Fig.HJshows that the Pt-dependence for 7r + and 7r~ are very similar to each 
other for each target, but that the slopes for the deuteron target are somewhat smaller than 
those for the proton. For a more quantitative understanding of the possible implications, 
we study the data in the context of a simple model in which the Pt dependence is described 
in terms of two Gaussian distributions for each case. Following Ref. we assume that the 
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FIG. 3: The parameters A and B [the relative coefficients of the cos<p (ctlt) and cos2</> (cttt) 
terms] averaged over tt + and ir~ detected from proton and deuteron targets, as a function of x 
at (z)=0.55 (left), and as a function of z at (x)=0.32 (right). The average value of transverse 
momentum ((|-Pt|)) is ~0.05 GeV. The dashed lines indicate the weighted averages for z < 0.7, 
which are also enumerated in each panel. Errors indicated include only statistical contributions. 
Systematic errors are highly correlated from point to point, and are estimated at 0.03 on both 
A and B. The open symbols are from exclusive pion production (see text). The solid lines are 
theoretical predictions 

widths of quark and fragmentation functions are Gaussian in k t and p t , respectively, and 
that the convolution of these distributions combines quadratically. The main difference from 
Ref. [6J is that we allow separate widths for up and down quarks, and separate widths for 
favored and unfavored fragmentation functions. The widths of the up and down distributions 
are given by \i u and fid, respectivel y, a nd the favored (unfavored) fragmentation widths are 
given by fi + (//_). Following Cahn and more recent studies [6(, we assume that only the 
fraction z of the quark transverse momentum contributes to the pion transverse momentum 
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FIG. 4: The -P t 2 dependence of differential cross-sections per nucleus for ir^ production on hydrogen 
(H) and deuterium (D) targets at (z}=0.55 and {x)=0.32. The solid lines show the result of the 
seven-parameter fit described in the text. The error bars are statistical only. Systematic errors 
are typically 4% (relative, see text for details). The average value of cos(<p) varies with P 2 (see 
Table 1. 

(see Fig. [1]). We assume further that sea quarks are negligible (typical global fits show less 
than 10% contributions at x — 0.3). To make the problem tractable, we take only the 
leading order terms in (kt/Q), which was shown to be a reasonable approximation for small 
to moderate P t in Ref. p. The simple model can then be written as: 

a ;+ = C[A Cl (P t )e^ p ? + (d/u)(D-/D + )c 2 (P t )e- b a p ?} 
a;~ = C[4(D-/D + )c 3 (P t )e- b " p ? + {d/u)c A {P t )e- b * p ?\ 
<+ = C[A(d/u) Ci (P t )e~ b t p i + (D-/D+)c 3 (P t )e- b " p ?} 
<- = C[A(d/u){D-/D+)c 2 (P t )e- b * p ? + c x {P t )e- b ^\ 

where C is an arbitrary normalization factor, and the inverse of the total widths for each 
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combination of quark flavor and fragmentation function are given by 



(5) 

6j = (^2 + 4)- 1 



and we assume = (a p + a n )/2. The Cahn effect {(], [ll| is taken into account through the 
terms: 

c 1 (P t ) = l. + c b (P t ,<cos(0) >){J*b+ 
c 2 (P t ) = 1. + Cb(P t , < cos(0) >)A*d&d 

c 3 (P t ) = 1. + Co(P t , < cos(0) >)fi 2 u b- (6) 
c 4 (P t ) = 1. + Co(P t , < cos(0) >)^6+ 



Co (p, < cos(0) >) = ^.y^y j? < cos(0) > . 

We fit for the four widths (fi u , fid, and C, and the ratios D~ / D + and d/tt, where 
the fragmentation ratio is understood to represent the data-averaged value at z — 0.55, 
and the quark distribution ratio is understood to represent the average value at x — 0.3. 
The fit describes the data reasonably well (x 2 = 78 for 73 degrees of freedom), and finds 
the somewhat low ratio d/u = 0.30 ± 0.03 (the LO GRV98 fit has about 0.40 for 
valence quarks), and the more reasonable ratio D~/D + = 0.42 ± 0.01 (a fit to HERMES 



results [22] , D~/D + = 1/(1 + z) 2 , predicts 0.42 at z = 0.55). Both d/u and D~ / ' D + are 
largely uncorrelated with other fit parameters. Since the data are at fixed z, the main 
contributions that distinguishes large fragmentation widths from large quark widths are the 
^-dependence Calm-effect q terms. While there is a significant inverse correlation between 
the two most important quark and fragmentation widths, (fi u and respectively), the fit 
finds a clear preference for \x u to be smaller than /x + as shown in Fig. [5K On the other hand, 
the fit finds fid and /i_ to be of the same magnitude and not strongly correlated, as shown 
in Fig. [5b. 

The fit tends to favor a larger k t width for d quarks (fi 2 , = 0.22 ± 0.13 GeV 2 ) than for u 
quarks (fi 2 = — 0.01±0.04 GeV 2 ), as illustrated in Fig. [5b, although the error on the d quark 



width is relatively large. The tendency is consistent with a di-quark model [23j in which the 
d quarks are only found in an axial di-quark, while the u quarks are predominantly found 
in a scalar di-quark. If the axial and scalar di-quarks have different masses, for example 
0.9 and 0.6 GeV, then the d quark distribution falls off more slowly with k t than the u 
quark distribution. In this model, the distributions show considerable deviation from an 
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FIG. 5: Fit parameters (crosses) and two-standard-deviation contours from the seven-parameter 
fit to the data shown in Fig. 3: a) u quark width squared // 2 versus favored fragmentation width 
squared fi^_; b) fi d versus u; jj c) // 2 versus d) /i 2 vs The large dot near the bottom of panel 



c is from a di-quark model 
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231 ] . The dashed line in panels c and d indicate ^ = ^ and ^ 2 = fJ,\, 



exponential falloff, but if we take the slope at the origin, the corresponding widths are 
/i 2 = 0.04 and = 0.08 GeV 2 , as illustrated by the solid dot in Fig. [5^. Fixing the quark 
widths to these values still gives a reasonable fit to our data (x 2 = 81 for 75 degrees of 
freedom). The magnitude of both widths is moderately sensitive to the choice of the model 
parameter Ao (we used 0.6 GeV), although the difference in widths is largely driven by the 
difference in axial and scalar di-quark masses (for example, increasing the axial di-quark 
mass to 1.4 GeV increases /i 2 , to 0.25 GeV 2 , which is the central value of our fit). Using the 
fit parameters, we find the magnitude of the cos(0) term A (Eq. [2]) at P t = 0.4 GeV to be 
about -0.05 for all cases except tt + production from hydrogen, where it is 0.01. These results 
are similar in sign and magnitude to those found in the HERMES experiment [24 ] . 

We find that the fragmentation widths /i+ and /i_ are correlated, as illustrated in Fig. [5]i, 
although the allowed range is not large, and the central values (fi 2 ^ = 0.23 ± 0.04 GeV 2 and 
/i 2 = 0.19 ± 0.04 GeV 2 ) are in reasonable agreement with both each other and also the 
flavor-averaged value of 0.20 GeV 2 found in Ref. 6j. While there is a slight tendency for the 
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favored width to be larger than the unfavored one, a reasonable fit can be obtained setting 
the widths equal to each other ( X 2 = 81 for 74 d.f., /i 2 + = /i 2 _ = 0.22 ± 0.04 GeV 2 ). 

To estimate the effect of experimental systematic errors on our fit results, we repeated 
the fits with: no diffractive p subtraction; no exclusive radiative tail subtraction; relative 
target thickness changed by 1%; and difference in 7r + and 7r~ absorption changed by 1%. In 
all cases, the quark and fragmentation width results remained well within the error ellipses 
shown in Fig. [51 The only parameter that changed significantly is the d/u ratio, which goes 
up to 0.33 with no exclusive tail subtraction. We found no significant change to the fit 
parameters upon adding to p 2 and /i 2 , an average nucleon transverse momentum squared of 
0.001 GeV 2 (evaluated using the Paris wave function [251 ]) for the deuteron model. 

In summary, we have measured semi-inclusive electroproduction of charged pions (jr ) 
from both proton and deuteron targets, using 5.5 GeV energy electrons at Jefferson Lab. We 
find the azimuthal dependence to be small, compared to exclusive pion electroproduction, 
and consistent with theoretical expectations • In the context of a simple model with 

only valence quarks and only two fragmentation functions, we find the k t width of d quarks to 
be larger than for u quarks, for which the width is consistent with zero within the statistical 
error. We find that the favored fragmentation p t width to be (unexpectedly) somewhat 
larger than the unfavored width, although both are larger than the u quark width. 

All of the above fit results can only be considered as suggestive at best, due to the 
limited kinematic range covered, the somewhat low u/d ratio that we find, and the very 
simple model assumptions described above. Many of these limitations could be removed 
with future experiments covering a wide range of Q 2 (to resolve additional higher twist 
contributions), full coverage in 0, a larger range of P t , a wide range in z (to distinguish 
quark width terms, which is weighted by powers of z, from fragmentation widths, which 
likely vary slowly with z), and including the ir° final state for an additional consistency 
check (particularly on the assumption that only two fragmentation functions are needed for 
charged pions from valence quarks). Some of these goals should be attainable with existing 
and upcoming data from Jefferson Lab, especially after the planned energy upgrade. These 
data should provide potential information on how hadron transverse momentum in SIDIS is 
split between fragmentation and intrinsic quark contributions. 
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TABLE I: Differential cross-sections per nucleus for ^ production on hydrogen and deuterium 
versus Pf. The error bars are statistical only. The values of cos(</>) and cos(20) averaged over the 
experimental acceptance are also indicated (see Fig. [2|). 

Pf < cos(c/>) > < cos(20) > of of of of 



GeV 2 nb/sr/GeV 3 nb/sr/GeV 3 nb/sr/GeV 3 nb/sr/GeV 3 
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